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Abstract 

Accurately quantifying soil organic carbon (SOC) Is considered fundamental to studying soil quality, modeling the global 
carbon cycle, and assessing global climate change. This study evaluated the uncertainties caused by up-scaling of soil 
properties from the county scale to the provincial scale and from lower-level classification of Soli Species to Soil Group, 
using four methods: the mean, median. Soil Profile Statistics (SPS), and pedologlcal professional knowledge based (PKB) 
methods. For the SPS method, SOC stock is calculated at the county scale by multiplying the mean SOC density value of 
each soil type in a county by Its corresponding area. For the mean or median method, SOC density value of each soil type is 
calculated using provincial arithmetic mean or median. For the PKB method, SOC density value of each soil type is 
calculated at the county scale considering soil parent materials and spatial locations of all soil profiles. A newly constructed 
1:50,000 soil survey geographic database of Zhejiang Province, China, was used for evaluation. Results indicated that with 
soli classification levels up-scaling from Soli Species to Soil Group, the variation of estimated SOC stocks among different soil 
classification levels was obviously lower than that among different methods. The difference in the estimated SOC stocks 
among the four methods was lowest at the Soil Species level. The differences In SOC stocks among the mean, median, and 
PKB methods for different Soil Groups resulted from the differences in the procedure of aggregating soil profile properties 
to represent the attributes of one soil type. Compared with the other three estimation methods (i.e., the SPS, mean and 
median methods), the PKB method holds significant promise for characterizing spatial differences In SOC distribution 
because spatial locations of all soil profiles are considered during the aggregation procedure. 
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Introduction 

Soil organic carbon (SOC), which plays a critical role in the 
global carbon cycle, comprises a major part of the terrestrial 
carbon reservoir [1-3]. In terrestrial ecosystems, SOC stock is 
almost three times the size of carbon storage in the vegetation of 
terrestrial ecosystems [4], and approximately twice as large as 
carbon storage in the atmosphere [5]. Because of the important 
role of SOC and its large quantity stored in terrestrial ecosystems, 
a slight change in SOC stock may influence global climate [6-8]. 
Accurately quantifying SOC stock has become a focus of present 
research on global climate change, and is considered essential for 
studying soil quality, modeling the global carbon cycle, and 
assessing global climate change [8-12]. 

Different SOC stock estimates, however, can vary gready at 
both global and regional scales [13-17]. For global scales, Bohn 
[18] estimated the total SOC stock was 3,000 Pg (1 Pg= lO''^ g), 
whereas Bolin [19] estimated only 710 Pg, over a four-fold 
difference. In China, estimated SOC stocks for terrestrial 
ecosystems range from 50 Pg [20] to 185.7 Pg [21], also 
approximately a four-fold difference. The methodology for 
estimating SOC stocks, such as soil profile-based [22-24] and 



model-based (e.g., the CENTURY [25] and DeNitrification- 
DeComposition models [26]), is one of the main factors attributed 
to the wide range of differences in SOC stock estimation from 
different studies. The soil profile-based methodology calculates 
SOC stock using soU profiles and their corresponding areas 
obtained from soil survey products (the soil type method), 
vegetation type maps (the vegetation type method), or life zone 
maps (the life zone method), among which the soil type method is 
the most widely used [9, 1 8- 1 9,27-28] . According to the sources of 
soil area, the soil type method contains the Soil Profile Statistics 
(SPS) method and the GIS-based Soil Type (GST) method [2]. 
The SPS method calculates SOC stock by multiplying the SOC 
density value of a soil type by its corresponding field survey area 
recorded in soU survey reports (e.g.. Soil Species of China [29]). 
The GIS-based Soil Type method calculates areas of various soil 
types accurately based on digital soil map and can provide 
information on the spatial distribution of SOC stocks. 

In China, different scale (i.e., county, provincial, and national) 
soil survey products (e.g., soil survey reports, soU maps) of the 
Second National Soil Survey of China are the most important data 
sources for SOC stock estimations [15,20-21,30-33]. For the 
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Figure 1. Study area and distribution of soil profile sites. 
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Second National Soil Survey of China, soil profile properties were 
aggregated sequentially in the county, provincial, and national 
scales. Different number of sampling profiles can be used to 
represent the attributes of one soil type at diflFerent scales and can 
cause uncertainties in SOC stock estimation. At the provincial 
scale, calculating the mean and median of SOC density values of 
multiple soil profiles of the same soil type name to represent the 
SOC density value of that soil type are the two most commonly 
used GIS-based Soil Type methods [32-34]. At the county scale, 
the SPS method [20,30-32] and the pedological professional 
knowledge based (PKB) method [2,15,34-36] are the two most 
commonly used methods. Compared with the median or mean 
method, the PKB method aggregates soil profile properties 
downscaling from the provincial scale to more detail soil map 
units, which links a SOC density value of each soil profile to a 
digital soil map according to the identity or similarity in soil parent 



materials and spatial locations of all soil profiles at the county 
scale. However, the uncertainties of soil profile properties caused 
by up-scaling from the county scale to the provincial scale among 
the four methods remain unknown. 

Soil map up-scaling includes soil classification level up-scaling 
and resolution up-scaling. The former aggregates Soil Species to 
higher soil classification levels (e.g.. Soil Group, Subgroup, Soil 
Family); the latter aggregates higher resolution soil maps to lower 
resolution ones, such as from 1 :50,000-county-scale map to 
1 : 1 ,000,000-country-scale map. Both the soil classification level 
up-scaling and resolution up-scaling can cause uncertainties in 
SOC stock estimation. Zhao et al. [34] tested the effects of soil 
map resolution up-scaling from 1:500,000 to 1:10,000,000 using 
the three GIS-based Soil Type methods (i.e., the mean, median, 
and PKB methods) on SOC stock estimation for Hebei Province, 
China. However, the effects of these three GIS-based Soil Type 
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Table 1. Methods used to estimate SOC (soil organic carbon) stocks in Zhejiang Province, China. 





SOC density 
Method value 


Scale 


Note 



SPS" Mean 



Mean or Mean or median 
median 



PKB" Mean 



County; soil (1) One soil species has multiple areas, which were surveyed county 

species by county. (2) Mean SOC density value was calculated from one 

or multiple profiles within the county. (3) 2154 soil profiles were used. 

Province; soil (1) One soil species has one area calculated from the digital soil map. 
species (2) Mean or median SOC density value was calculated from one 

or multiple profiles within the Province. (3} 2154 soil profiles were used. 

County; soil (1) Soil map units were derived from the digital soil map county 

map unit by county. (2} One soil map unit in one county may have one 

or multiple areas calculated from the digital soil map. (3) 
Mean SOC density value was calculated from one or 
multiple profiles located within one polygon; polygons 
belong to one soil map unit in one county may 
assigned different SOC density values. (4) 
2154 soil profiles were used. 



^SPS, Soil Profile Statistics; PKB, pedological professional knowledge based. 
doi:1 0.1 371 /journal.pone.0097757.t001 



methods on SOC stock estimation at soil map resolution larger 
than 1:500,000 remain unknown. Moreover, the uncertainties in 
SOC stock estimation caused by soil classification level up-scaling 
using the SPS, mean, median, and PKB methods remain 
unknown. 

The goal of this study is to evaluate the uncertainties of changes 
in scale among the four estimation methods (i.e., the SPS, mean, 
median, and PKB methods) by using a newly completed 1 :50,000 
soil survey geographic database of Zhejiang Province, China. 
Specifically, the main objectives of this study were to: (1) evaluate 
the uncertainties caused by up-scaling of soil profile properties 
from the county scale to the provincial scale and soil classification 
levels from Soil Species to Soil Group among the four methods; 
and (2) quantify spatial differences in SOC stock estimation among 
the three GIS-based Soil Type methods. 



km^. Hills and mountains are dominant terrain, occupying 70.4% 
of the total land area. The plains and basins account for 23.2%, 
and the remaining 6.4% is comprised of lakes and rivers. 
Provincial topography is characterized by high mountains in the 
southwest and low plains in the northeast. The elevation ranges 
from 0 to 1895 m with an average of 296 m above sea level. 

With a subtropical monsoon climate, Zhejiang Province has an 
annual average temperature from 15 to 18°C and an annual 
precipitation of 1200 to 1800 mm. As substantial differences exist 
in climate, geomorphology, geohydrology, land use, and parent 
material throughout the Province, its vegetation and soils and their 
spatial distribution patterns vary gready. According to Soils of 
Zhejiang Province [37], Red soils are the dominant Soil Group, 
accounting for approximately 39% of the total area of Zhejiang 
soils. 



Materials and Methods 

Study Area 

Zhejiang Province is located between 27°06' and 31°03' N 
latitude and 118°01' and 123°10' E longitude in southeastern 
China (Figure 1). The Province covers a land area of 1.02*10'^ 



1100 




Soil Group Subgroup Soil Family Soil Species 

Figure 2. The estimates of SOC (soil organic carbon) stocks for 
soil classification levels up-scaling from Soil Species to Soil 
Group using the mean, median, SPS (soil profile statistic), and 
PKB (pedological professional knowledge based) methods. 

doi:1 0.1 371 /journal.pone.0097757.g002 



Data Sources 

A recently completed 1:50,000 soil survey geographic database 
of Zhejiang Province [38] was used in the study, which is the most 
detailed soil survey data at the provincial scale in China to date. 
The database includes Soil Spatial Data and Soil Attributes Data. 
The Soil Spatial Data include a 1 :50,000 digital soil map of the 
Province and a digital map of 2 1 54 geo-referenced soil sampling 
profile sites, both of which were derived by digitizing and re- 
compiling field soil survey maps at 1:50,000 from 76 counties in 
Zhejiang. The Soil Attributes Data contain the properties of the 
2154 geo-referenced soil sampling profiles, which were taken from 
original county soil survey reports. For each profile, there are one 
to seven soil layers; for each layer, there are up to 1 04 soil physical 
and chemical properties, including geographic location, depth, 
bulk density, organic matter content, and gravel content, etc. 

The soil survey geographic database was used to estimate SOC 
stocks with the mean, median, and PKB methods. This 
information from the Second National Soil Survey of China 
conducted in early 1980's was the most comprehensive and 
detailed inventory of soil characteristics in Zhejiang Province. The 
Genetic Soil Classification System of Zhejiang Province was used 
during the field soil survey. The soils were classified using a six- 
level hierarchical scheme: Order, Suborder, Group, Subgroup, 
Family, and Species. The Soil Species is the basic classification 
level, and the Soil Group is the most commonly used classification 
level in China [2]. According to the Genetic Soil Classification 
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System of Zhejiang Province, the collected 2154 soil profiles 
originated from 10 Soil Groups, 21 Subgroups, 99 Soil Families 
and 277 Soil Species. 

Estimation of Organic Carbon Density for Soil Profiles 

For each of the 2 154 soil profiles, SOC density was calculated 
with the formula [34,39]: 



I SOCDz,= ^[(l-6,%)xp,.xC,xr,/100] (1) 

S ' = 1 

c 

^ where SOCD^ (kg m ^) is the SOC density of a soU profile within 

c a depth of D (cm), n is the number of soil layers in the soil survey, 

S 6,% represents the volumetric percentage of gravel (>2 mm) 

> content, p; is the soil bulk density (g/cm"'*), C, is the organic carbon 

° content (C g/kg), and TJ represents the thickness (cm) of the layer i. 

S Density of SOC was estimated to a maximum depth of 100 cm to 

it facilitate comparison among data sets. For profiles with actual 

S depths greater than or equal to 1 00 cm, but less than 1 00 cm was 

.1 observed, data for the unobserved profile sections were derived 
from the mean values of all the corresponding soil profiles of the 

. same Soil Species or Soil Family [15,33,40]. Organic carbon 

E S content is calculated by multiplying soil organic matter content by 

o S 0.58 (the Bemmelen index), which is based on the assumption that 

1^ g- soil organic matter contains approximately 58% organic carbon 

3 i m- 

f i 



"S E 



Estimation of Regional Soil Organic Carbon Stocks 

The SOC stocks of the Province were estimated by using the 

2 I SPS, mean, median, and PKB methods (Table 1). 

j£ For the SPS method, the mean SOC density value of each soil 

<g 2 type in one county was calculated from one or multiple profiles 

S = belonging to that soil type; then the SOC stock value for that soil 

I type was calculated by multiplying its SOC density value by its 

.£ corresponding area; finally, the SOC stock values of all soil types 

1 .1" were summed up as the total SOC stock for the Province. The 
^ g area calculated from the 1 :50,000 digital soil map was used in this 

° Study to facilitate comparison among dilferent methods. 

E ■= For the mean or median method, the arithmetic mean or 

Si o . . . 

>, Q. median of SOC density values of all soil profiles of the same soU 

g o type name in the Province was calculated fu-st as the SOC density 

2 o value of that soil type; secondly, SOC stocks for polygons on the 
^ .Q 1 :50,000 digital soil map were calculated by multiplying the mean 
(5 I or median values by corresponding areas of the soil types and then 
^ <u were summed up as the total stock for the Province. 

'5 For the PKB method, the SOC density value of each of the 

2154 soil profiles was linked to corresponding polygons on the 

1:50,000 digital soil map according to the same soil type name, the 

^ J! soil parent materials, and the spatial locations by county. When 

^ -g two or more soil profiles of the same soil type name were located in 

£ J ^ one polygon, the mean SOC density value of these soil profiles was 

= !; 8 calculated and used for the linkage. The SOC stock values for all 

"2 e ^ polygons were calculated and then were summed up as the total 

I ^ S; stock for the Province. 
5 g o 

^ _^ ° The estimation of SOC stocks using the three GIS-based Soil 

g ^- ° Type methods were performed using ESRI software ArcGIS 10.0 

5 ^. 1 (Redlands, CA) with water and urban areas excluded from the 

S 5 calculation. 
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SOC stock (10' t) 
■ 0-20 cm g 20-40 cm ^ 40-60 cm 0 60-80 cm g 80-1 00 cm 

Figure 3. Estimated SOC stocks at different soil depths using 
the mean, median, SPS (soil profile statistic), and PKB 
(pedological professional knowledge based) methods. 

doi:1 0.1 371 /journal.pone.0097757.g003 

Results 

Effect of Soil Classification Levels Up-Scaling from Soil 
Species to Soil Group on SOC Estimation Using Different 
Methods 

The SOC stocks for soil classification levels up-scaling from Soil 
Species to Soil Group using the mean, median, SPS, and PKB 
methods were calculated (Figure 2). Among the four methods, the 
estimated SOC stocks for the study area using the median method 
were always the lowest regardless of soil classification level. 
Because the SOC densities of the 2154 soil profiles showed a 
positively skewed distribution with skewness of 1 4. 1 7 (Table 2), the 
estimated SOC stocks using the mean method were much higher 
than those using the median method at all four soil classification 
levels. 

With soil classification levels up-scaling from Soil Species to Soil 
Group, the estimated SOC stocks presented small variations with 
coeflTicient of variations of 2.1%, 1.4%, 1.4%, and 1.2% for the 
mean, median, SPS, and PKB methods, respectively. The 



coefficient of variations of estimated SOC stocks among different 
methods presented an increasing trend, with 4.4%, 6.2%, 6.3%, 
and 7.1% for soil classification levels of Soil Species, Soil Family, 
Subgroup, and Soil Group, respectively. The estimated SOC 
stocks using the mean, median, and SPS methods at the SoU 
Species level were 3.1%, 10.6%, and 0.5% lower than that 
obtained by the PKB method, respectively. The variation of 
estimated SOC stocks among different methods was obviously 
larger than that among different soil classification levels. 

Effect of Different Methods on the Estimates of SOC 
Stocks for Various Soil Groups 

Soil Group is the most stable and consistent soil classification 
level commonly used in China [2]. In this study, the SOC stocks 
and SOC densities for various Soil Groups estimated by the mean, 
median, SPS, and PKB methods are presented in Table 3. For the 
SoU Groups with SOC density values of SoU profiles presented 
positively skewed distribution, the estimated SOC stocks using the 
median method were usually lower than that using the mean and 
PKB methods (Table 2, Table 3). The estimated SOC stocks for 
various Soil Groups using the SPS method were most similar to 
that obtained by the PKB method because both methods 
aggregate soil profile properties at the county scale. The estimated 
SOC stocks for various SoU Groups using the mean and median 
methods showed large differences compared with the PKB method 
because both the mean and median methods aggregating soil 
profile properties at the provincial scale that the SOC density 
values of all soil profiles belonging to a soil type in the Province 
were aggregated to one value. 

The most pronounced differences of the estimated SOC stocks 
for the ten Soil Groups (i.e.. Red soils. Yellow soils. Purple soils, 
Limestone soUs, Skel soils. Red clay soils. Mountain meadow soils, 
Fluvio-aquic soUs, Coastal saline soUs and Paddy soUs) among the 
three GIS-based SoU Type methods occurred on YeUow soUs and 
Mountain meadow soils, which presented high SOC density values 
with high coefficient of variances. Thus, the process of linking the 



A 




150 200 



Mean vs. PKB 
Difference of SOC 
density (kg m-2) 

< -50 -2-2 
-50 --10 2-5 
-10 --5 ^^5- 10 
-5- -2 ^■>10 

Water and urban area 




0 25 50 



100 



Median vs. PKB 
Difference of SOC 
density (l<g m-2) 

< -50 -2-2 
-50 --10 2 -5 

-10--5 ^^5-10 
-5- -2 ^■>10 

Water and urban area 



Figure 4. Spatial differences of estimated SOC (soil organic carbon) densities between the mean and PKB (pedological professional 
knowledge based) methods (A), and between the median and PKB methods (B). 

doi:1 0.1 371 /journal.pone.0097757.g004 
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aggregated SOC density values of soil profiles belonging to these 
two Soil Groups to corresponding polygons on the 1 :50,000 digital 
soil map would significandy influence the estimates of SOC stocks. 
One important factor leading to the significantly high variation of 
SOC in Mountain meadow soils was the limited number of soU 
profiles (n = 4) in the soil survey area. 

Effect of Different Methods on the SOC Stock Estinnation 
at Different Soil Depths 

For the four methods, the trend of estimated SOC stocks at each 
of the- fi\-c soil depths (0-20, 20-40, 40-60, 60-80, and 80- 
100 cm) is consistent with the soil depth of 0-100 cm, followed the 
order: PKB>SPS>mean>median (Figure 3). At each of the five 
soil depths, the estimated SOC stock using the SPS method 
presented the smallest difference while the median method 
presented the largest difiFerence compared with the PKB method. 
Decreasing SOC stock with increasing soil depth was evident 
independent of estimation method, and these decreasing trends 
are consistent among aU four methods. At the depth of 0-20 cm, 
SOC stocks estimated by the mean, median, SPS, and PKB 
methods account for 41.5%, 42.3%, 41.1%, and 41.1%,, of total 
SOC stocks, respectively, while the SOC stocks at the depth of 80— 
100 cm estimated by the four methods account for only 9.2%, 
8.9%, 9.1%, and 9.1%, of die total SOC stocks, respectively. 

Comparison of the Differences in the Spatial Patterns of 
SOC Distribution among the Three GIS-based Soil Type 
Methods 

Estimated SOC densities using three GIS-based Soil Type 
methods showed clearly spatial differences (Figure 4). The largest 
differences were as high as —229.63 kg m ~ between the median 
and PKB methods and —175.20 kg m between the mean and 
PKB methods. Polygons mosdy belonging to the Soil Groups of 
Mountain meadow soils and Yellow soils showed substantial 
differences (< — 50 kg m~^) in SOC densities among the three 
GIS-based Soil Type methods. For the polygons with large spatial 
difiFerences in SOC densities ranging from —50 kg m ^ to —5 kg 
m~^ and >5 kg m~^, the locations were usually in mountain and 
hilly areas with dominate Soil Groups of Yellow soils, Red soils, 
and Skew soils, with the remainder in plain areas with dominate 
Soil Group of Paddy soils. The total areas for spatial dififerences in 
SOC densities with absolute value larger than 5 kg m~^ between 
the mean and PKB methods and the median and PKB methods 
were estimated to be 12883.4 km^ and 11428.0 km^, accounting 
for 12.8% and 11.3% of the total area of Zhejiang soils, 
respectively. The total areas for spatial difiFerences in SOC 
densities with absolute value lower than 2 kg m~^ between the 
mean and PKB methods and the median and PKB methods were 
estimated to be 55316.1 km^ and 59147.7 km^, accounting for 
54.9%i and 58.7% of the total area of Zhejiang soils, respectively. 

Discussion 

With soil classification levels up-scaling from Soil Species to Soil 
Group, the estimated SOC stock values using the SPS method 
were most similar to that obtained by the PKB method. The 

difference of the estimated SOC stocks between the SPS method 
and the PKB method occurred in the cases where two or more soil 
profiles belonging to one soil species existed in one county. 
However, for most cases, there was only one soil profile belonging 



to one soil species in one county, the SOC density value of that soil 
profile was used to represent the SOC density value of that soil 
species in that county, which resulted in the same estimated SOC 
stock value of that soil species in that county using the SPS and 
PKB methods. Therefore, more number of soil profiles is needed 
to further compare the difference between these two methods. 

The estimated SOC stock values using the mean method were 
similar to that using the PKB method. When two or more soil 
profiles located in the same polygon within a county, the mean 
SOC density value of these soil profiles was calculated and used to 
link with that polygon, so the PKB method approaches the mean 
method in Unking SOC density values with polygons to some 
extent [34], especially for higher soil classification levels (e.g.. Soil 
Group, Subgroup, Soil Family). 

With soil profile properties up-scaling from the county scale (the 
SPS and PKB methods) to the provincial scale (the mean and 
median methods), the SOC density values of all soil profiles 
belonging to a soil type in the Province were aggregated to one 
value. The aggregating procedure for the mean or median method 
could cause large uncertainties because information about how 
much area one profile can represent was missing, thus difiFerent 
polygons belonging to the same map unit on the digital soil map 
were linked with the same SOC density value. For the PKB 
method, different polygons belonging to the same map unit may 
be assigned difiFerent SOC density values depending on their 
locations, thus has an obvious advantage in the demonstration of 
spatial differences in SOC distribution. The difiFerence in 
aggregating procedure led to Soil Groups (e.g., Mountain meadow 
soils. Yellow soils) with soil profiles of high variations in SOC 
density values presented significant spatial difiFerences among the 
three GIS-based Soil Type methods. 

Conclusions 

The up-scaling of soil classification levels from Soil Species to 
Sou Group has small effect on SOC stock estimation. However, 
obvious differences occurred with difiFerent estimation methods, 
especially for Soil Groups with high variations in SOC densities 
and spatial distribution. The PKB method, which links soil profile 
properties to spatial databases by county, produced more stable 
results for soil types, soil parent materials, and spatial locations of 
aU soil profiles under consideration. Thus this method has an 
obvious advantage in the demonstration of differences in spatial 
patterns of SOC distribution than both the mean and median 
methods. We recommend the PKB mc-thod as a prior option 
rather than the mean, median, and SPS methods for SOC stock 
estimation in China, especially when 1:50,000 soil survey 
geographic database (soil map) is available. It potentially reduces 
uncertainties related to up-scaling soil profile properties. 
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